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Scattering theory traditionally deals with the asymptotic behaviour of a system far removed from 
the actual scattering event. Here we present an experimental study of the one-dimensional scattering 
of a non-interacting condensate of 87 Rb atoms from a potential barrier in the non-asymptotic regime, 
for which the collision dynamics are still ongoing. We show that for near-transparent barriers, 
there are two distinct transient scattering effects that arise and dramatically change the momentum 
distribution during the collision: a deceleration of wavepacket components in mid-collision, and an 
interference between incident and transmitted portions of the wavepacket. Both effects lead to the 
re-distribution of momenta giving rise to a rich interference pattern that can be used to perform 
reconstruction of the single-particle phase profile. 
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Scattering is traditionally described asymptotically, 
connecting the initial and final scattering states, while 
bypassing the dynamics during the collision itself [1]. 
In typical particle scattering experiments, measurements 
are made long after the collision and thus the asymptotic 
solutions provide an accurate description of the observa- 
tions. Yet, there are circumstances in which the scatter- 
ing cannot be described asymptotically and knowledge of 
the full wavefunction is required, for example in under- 
standing the quantum kinetics of moderately dense gases 
[2, 3]. Brouard and Muga [4, 5] theoretically studied the 
one-dimensional scattering of a wavepacket from a delta- 
function potential. They showed that during the collision 
with a repulsive potential, the momentum-space wave- 
function could exhibit dramatically different features, for 
instance the generation of high- momentum components. 
This non-classical momentum enhancement is a conse- 
quence of the wave nature of matter, in the spirit of 
quantum reflection [6-8] - an effect that occurs when 
the potential changes abruptly on the scale of the de 
Broglie wavelength, irrespective of the sign of the change. 
However, the effect described here is distinctly transient, 
manifesting itself only during the scattering event, and 
vanishing in the asymptotic limits. 

Over the past few decades, impressive advances in ex- 
perimental techniques have granted access to ever faster 
timescales and lower energies, allowing for direct probes 
of the no n- asymptotic scattering regime. For example, 
ultrafast laser pulses can now be used to probe sub- 
femtosecond timescales, providing time-resolved probes 
of electron dynamics [9]. In parallel, atom cooling tech- 
niques [10-12] have matured, and now routinely produce 
samples in the nano-Kelvin regime [13, 14]. The scat- 
tering dynamics of these ultracold systems occur on an 
easily accessible microsecond timescale, making them an 
ideal system for investigation of non-asymptotic scat- 
tering. Bose-Einstein condensates of dilute gases have 



been used extensively to observe matter- wave phenomena 
[8, 15, 16]. Due to the relatively low densities, the inter- 
particle interactions are weak, allowing for clear observa- 
tion of single-particle quantum effects. Furthermore, the 
dynamic Stark effect allows for creation of nearly arbi- 
trary potentials with spatial features limited only by the 
wavelength of light used to induce them [17], and can 
be modulated much faster than the dynamics of the con- 
densate [18]. Here we make use of these advantages to 
experimentally study scattering of a matter wavepacket 
from an optically-induced potential barrier in the non- 
asymptotic regime. We observe the momentum-space 
wavefunction reshaping effect predicted by Brouard and 
Muga, and further show that scattering from a finite- 
sized potential gives rise to a second effect due to the 
deceleration of wavepacket components in mid-collision. 
These transient scattering effects are observed to result 
in a rich momentum-space interference pattern. 

At first glance, the predicted momentum redistribu- 
tion and enhancement of high momenta might be inter- 
preted as a spatial compression of the wavepacket as it 
impinges upon the barrier, resulting in a broadened mo- 
mentum distribution. Indeed, it has been shown that 
the enhancement increases with the height of the bar- 
rier, saturating in the opaque barrier limit [5]. However 
Perez, Brouard and Muga [19] later showed that a much 
more significant momentum enhancement unexpectedly 
occurs for wavepackets with kinetic energy well above the 
barrier height. They interpreted this effect as an inter- 
ference between incident and transmitted portions of the 
wavepacket. In the limit of a nearly transparent barrier, 
the primary effect on the wavepacket is to write a phase 
shift. During the collision, half of the wavepacket has 
yet to reach the barrier and half has been transmitted. 
Thus there exists a phase discontinuity dividing the in- 
cident and transmitted components of the wavepacket in 
mid-collision (Fig. 1 inset). This sharp spatial feature in 
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FIG. 1. Transient enhancement of momentum during scatter- 
ing with a delta-function barrier. Shown is the simulated mo- 
mentum distribution, in the wavepacket centre-of-mass frame, 
long before/after the collision (red dashed line) and during 
the collision (blue solid line). Inset: Position- space distribu- 
tion during scattering (blue solid line) , and phase-profile (red 
dashed line). The barrier height has been chosen such that it 
results in a 7r-phase shift across the wavepacket. 



the phase results in destructive interference of the cen- 
tral momentum component and constructive interference 
in the tails [20] . The net result is a symmetrically broad- 
ened momentum distribution which exhibits non-classical 
momentum enhancement, and tails that fall off as 1/k 2 
(Fig. 1). For an incident wavepacket with linear chirp, 
as would result from free expansion, the momentum dis- 
tribution during scattering exhibits a much richer struc- 
ture. In mid-collision, the phase discontinuity is located 
at the spatial centre of the wavepacket. Only the central 
momentum components, predominantly located near the 
centre of the spatial wavepacket, acquire the long momen- 
tum tails. The high and low momentum components, 
which are spatially far from the discontinuity, remain 
unaffected. The resulting momentum distribution dis- 
plays rapid fringes due to interference between the redis- 
tributed central momentum component, and components 
in the tails of the distribution (Fig. 2(a)). 

For finite- width barriers, the interference effect is sig- 
nificantly reduced. This reduction depends exponen- 
tially on the product of the barrier width a y and the 
wavepacket momentum width Afc. In addition to the 
reduction of the momentum enhancement, there is a sec- 
ondary effect which is that a centre-of-mass momentum 
is imparted. Components of the wavepacket in mid- 
collision receive a net deceleration. These components 
are effectively pushed to lower momenta, resulting in in- 
terference between the incident and pushed components 
of the wavepacket (Fig. 2). This breaks the symme- 
try of the momentum distribution and is distinct from 
the interference between the transmitted and incident 
components. For our typical experimental parameters, 
<j y • Ak « 4, the transmitted-incident (T-I) interference 
effect is strongly suppressed, and the scattering is dom- 
inated by pushed-incident (P-I) interference. However, 
since the latter only changes the distribution of slow- 
moving components of the wavepacket and the interfer- 
ence increases our sensitivity to small amplitudes, it is 
still possible to unambiguously observe the T-I interfer- 
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FIG. 2. Momentum- space interference during scattering for a 
linearly chirped wavepacket. Simulated momentum distribu- 
tions, in the centre-of-mass frame of the incident wavepacket, 
before/after the collision (red dashed line) and during the 
collision (blue solid line), (a) The left column shows the de- 
pendence on barrier width. For a delta-function barrier, T-I 
interference is dominant, creating a symmetric distribution. 
For a finite- width barrier, P-I interference becomes dominant, 
breaking the symmetry, (b) The right column shows the effect 
of increasing the barrier height. 



ence effect on the fast-moving, transmitted components. 
An additional consequence of a finite-width scattering 
potential is that the strict 27r-phase periodicity in T-I in- 
terference that would occur for a delta-function potential 
is washed out by the continuous range of phases written 
on the wavepacket. 

It is important to emphasize the transient nature of 
both T-I and P-I interference effects. For our implemen- 
tation, the barrier height is typically 8 times smaller than 
the average kinetic energy of the wavepacket; there is no 
classical reflection and a probability of quantum reflec- 
tion less than 10 -13 . Yet even in this near-transparent 
barrier limit, there is a dramatic modification of the mo- 
mentum distribution during interaction with the barrier, 
which afterwards vanishes, returning to the original dis- 
tribution. In contrast, the only classical effect would be 
a temporary decrease in the momentum of the fraction 
of atoms in mid-collision at any given moment. 

We prepare a Bose-Einstein condensate of 87 Rb atoms 
in the ground state of an optical dipole trap formed by the 
intersection of two focused 980 nm beams (red-detuned 
from the 780 nm D2 transition) oriented in the horizon- 
tal x-z plane. The trap is cylindrical, with approximate 
radial and axial trapping frequencies uo x = uo y = 2tt • 300 
Hz and uo z = 2tt • 100 Hz, respectively. In this trap, we 
prepare nearly pure condensates of about 10 5 atoms. Af- 
ter preparation, we abruptly turn off the optical trap, 
dropping the atoms under gravity (^-direction) onto an 
optically-induced barrier potential. The barrier is posi- 
tioned beneath our atoms such that the collision typi- 
cally occurs after t co \ ~ 4 ms of free-fall, during which 
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the wavepacket acquires a linear chirp. The chirped 
wavepacket then collides with the barrier. To study the 
momentum distribution during the collision, we abruptly 
turn the optical barrier off within 20 /is, freezing out the 
dynamics of the collision. We then perform a long time- 
of-flight (TOF) expansion, followed by absorption imag- 
ing after a total time Uot — 30 ms. This TOF maps the 
momentum of the particles during the collision to their 
final position, such that the imaged position distribution 
is representative of the momentum distribution of the 
atoms during the collision. 

The barrier potential is generated by a 780 nm beam 
blue-detuned from the 87 Rb D2 transition by 150 GHz. 
This beam is approximately gaussian in shape with a 
large aspect ratio, and propagates along the x-axis. The 
rms barrier width along the direction of wavepacket prop- 
agation is <j y = 1.1 /im. The beam has a size <j z ~ 10 /im 
transverse to the wavepacket propagation, and a Rayleigh 
range of xr ~ 20 /im. a z is comparable to the size of the 
expanding cloud at the time of collision. We broaden the 
potential along z by rapidly scanning the position of the 
barrier using an acouto-optic deflector (AOD), extending 
the length of the barrier to ±60 /im. This scan occurs 
at 100 kHz and is much faster than the motion of the 
atoms, resulting in a time-averaged potential [18]. We 
tailor the AOD scanning waveform to generate a time- 
averaged potential flat to within 1% of the barrier height, 
ensuring that the collision dynamics are essentially one- 
dimensional. The barrier height is typically around fc^xl 
/iK, where fc# is Boltzmann's constant, and is chosen 
such that the effect of the barrier is to write a typical 
phase shift of approximately 37r, while remaining small 
compared to the typical wavepacket centre-of-mass ki- 
netic energy of ks x 8 /iK acquired in free-fall. 

Note that for an interacting system, such a phase im- 
printing would result in soliton formation [21, 22]. The 
effect studied here is a single-particle interference effect. 
The repulsive inter-particle interactions in the system are 
only relevant during the first millisecond of the experi- 
ment, driving the self-similar expansion of the cloud and 
thus determining the momentum distribution prior to 
collision. The mean-field energy is rapidly converted to 
kinetic energy on a timescale of l/oj y [23], and is negli- 
gible during the collision itself. This has been confirmed 
by comparing an interacting 3D Gross-Pitaevskii equa- 
tion simulation to a non-interacting Schrodinger equation 
simulation of the collision with matched pre-collision mo- 
mentum distributions. 

We measure the condensate momentum distribution 
during scattering for a variety of scattering parameters. 
Figure 3 shows a representative set of absorption im- 
ages and the corresponding integrated ID profiles for 
pre-collision expansion time t co \ =3.3 ms and increas- 
ing barrier height. The barrier height is expressed in 
terms of the estimated phase shift imparted to the fully 
transmitted atoms. As the barrier height is increased we 
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FIG. 3. Transient momentum-space interference during a col- 
lision. Representative absorption images and corresponding 
ID density profiles of a condensate after collision with a bar- 
rier at tool = 3.3 ms and TOF ttot = 30 ms. The barrier 
height is expressed in terms of the phase shift imparted to 
fully transmitted atoms. The cloud is moving under the in- 
fluence of gravity in the +?/-direction. The images have been 
low-pass filtered in the transverse ^-direction for presentation. 



observe the development of a rich interference pattern, 
consistent with P-I interference. For higher barriers, the 
central momentum component is more strongly deceler- 
ated, thus increasing the range over which interference is 
observed. The fringe pattern that develops is reflective 
of the momentum-space phase profile of the condensate 
prior to collision. For condensate expansion, we expect 
a quadratic phase profile [23], and a fringe spacing that 
decreases linearly from cloud centre. Note that our imag- 
ing system has a resolution of 3 /im, reducing the fringe 
visibility near the edges of the distribution. 

Given the finite barrier width, T-I interference is 
strongly suppressed and the scattering is dominated by 
P-I interference. However the latter effect is asymmetric, 
only affecting the low-momentum side of the distribution; 
thus T-I interference can still be observed on the high- 
momentum side. In our parameter range, the effect is 
most visible for large phase curvatures, which is achieved 
with long pre-collision times. Figure 4 shows a represen- 
tative ID profile (averaged over 20 images) for t co \ = 6.3 
ms and an estimated barrier phase shift of 3.17T. The 
experimental profiles show excellent agreement with our 
simulation, which accounts for the finite imaging resolu- 
tion and camera pixel size. Some jitter in the shot-to-shot 
position of the cloud may contribute to the reduced fringe 
visibility in the averaged experimental profiles when com- 
pared to simulation. At this pre-collision time, we begin 
to see T-I interference fringes on the high-momentum 
side. In principle, this type of interference also results 
in momentum enhancement; in our parameter regime, 
however, this effect is very small. Though the T-I inter- 
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FIG. 4. Evidence for transmitted- incident momentum-space 
interference for t co i — 6.3 ms, ttot = 30 ms, and an esti- 
mated barrier phase shift of 3.1 7r. (a) Sample experimental 
ID profile (average of 20 images), and (b) simulation (includ- 
ing image resolution effects). Shown are the density profiles 
for un-collided (red) and collided clouds (blue), and the dif- 
ference between these profiles, A. The pushed-incident in- 
terference is visible on the low- momentum side (negative y), 
while the transmitted-incident interference is visible on the 
high-momentum side (positive y) around y = 50 /im. 



ference fringes are clearly observed, the expected momen- 
tum enhancement signal is too small to be quantitatively 
extracted from the experimental noise. 

The observed momentum distributions during the col- 
lision exhibit a rich structure, containing both ampli- 
tude and phase information about the single-particle 
momentum-space wavefunction prior to collision. Both 
T-I and P-I effects result in interference between mo- 
mentum components initially at k\ and scattered by the 
barrier at t co \ to fc/, and un-scattered components which 
have momentum kf. The interference depends on the 
relative phase accumulated along the respective trajec- 
tories. By studying the fringe pattern we can extract the 
phase profile of the single-particle wavefunction. To illus- 
trate, we estimate the phase profile of our expanding con- 
densate (Fig. 5) by marking the maxima and minima of 
each P-I interference feature in the observed profiles and 
assigning a relative 7r-phase shift between adjacent fea- 
tures. For condensate expansion, a quadratic phase pro- 
file with a curvature that increases with the pre-collision 
time is expected for both momentum- and position-space 
distributions [23]. The position-space phase profile after 

TOF is 4>{z) = az\ where a * § - ^) . We 

find close agreement between a and the phase curvature 
extracted from simulated data using the above technique. 

Our ability to focus on the atom cloud centre after 
TOF is limited to an accuracy of ±300 /im. Off-centre 
imaging results in distinct diffraction features, and a gen- 
eral reduction in extracted experimental fringe curvature 
from the expected value. These imaging artifacts are 
well understood and entirely reproduced by our imaging 
simulations (performed after data collection). Thus we 
compare our extracted fringe curvature to that predicted 
from our simulations for a range of imaging planes from 
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FIG. 5. Extracting the phase profile from experimental data. 

(a) Sample profile for t co i = 3.3 ms, where each peak/trough 
is marked and a relative tt phase is assigned, (b) Each such 
phase profile is fitted to a quadratic, (c) The fitted phase 
curvatures are plotted against pre-collision time. The thick 
black line indicates expected free-expansion curvature, a. Red 
lines indicate the curvature extracted from simulated data for 
imaging on cloud centre, and increasingly off-centre (top to 
bottom; see text). The colour and symbols of fringe data in 

(b) match that of the fitted curvatures in (c). 



on-centre to 300 um off-centre. Note that the t co \ = 6.3 
ms data set where the T-I interference effect is most vis- 
ible (Fig. 4), is consistent with imaging on cloud centre. 

The technique used here to reconstruct the expanding 
condensate phase profile can be used for measurement 
of arbitrary single-particle phase profiles. This inter fero- 
metric approach bears some similarities to existing con- 
densate phase profile reconstruction techniques [16, 24]. 
In contrast to these approaches, our technique in prin- 
ciple allows for single-shot reconstruction, provided the 
scattering potential is sharp enough. For sharp scatter- 
ing potentials, T-I interference dominates, symmetrically 
scattering the central momentum component across the 
cloud. Thus the central momentum component acts like 
a local oscillator, beating against the other momentum 
components of the condensate simultaneously. On the 
other hand, for broad scattering potentials the P-I inter- 
ference dominates. Assuming a symmetric phase profile, 
it is still possible to extract the full phase information, 
as we have shown for condensate expansion. When us- 
ing P-I interference, detailed knowledge of the scattering 
potential and its effect on the wavepacket is required to 
accurately recover the relative phase between momentum 
components. Lastly, we note that since the fringe visibil- 
ity is linked to the off-diagonal elements of the density 
matrix, this technique can be extended to perform full 
tomography on the single-particle state. This could al- 
low for the characterization of many body-correlations in 
ultracold systems, for example during the superfluid to 
Mott insulator transition. 
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